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Abstract 

It is shown that the next-to-leading order (NLO) corrections to the QCD Pomeron 
intercept obtained from the BFKL equation, when evaluated in non-Abelian physical 
renormalization schemes with BLM optimal scale setting do not exhibit the serious prob- 
lems encountered in the MS-scheme. A striking feature of the NLO BFKL Pomeron 
intercept in the BLM approach is its rather weak dependence on the virtuality of the 
reggeized gluon. This remarkable property yields an important approximate conformal 
invariance. The results obtained provide an opportunity for applications of NLO BFKL 
resummation to high-energy phenomenology 
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The discovery of rapidly increasing structure functions in deep inelastic scattering 
(DIS) at HERA []J at small-x is in agreement with the expectations of the QCD high- 
energy limit. The Balitsky-Fadin-Kuraev-Lipatov (BFKL) resummation of energy 
logarithms is anticipated to be an important tool for exploring this limit. The leading 
order (LO) BFKL calculations predict a steep rise of QCD cross sections. Namely, 
the highest eigenvalue, u max , of the BFKL equation is related to the intercept of 
the Pomeron which in turn governs the high-energy asymptotics of the cross sections: 
a ~ s a ip- 1 = s ^ max _ The BFKL Pomeron intercept in the LO turns out to be rather large: 
otjp — 1 = uj™ ax = 12 In 2 (as/ft) — 0.55 for as = 0.2; hence, it is very important to know 
the next-to-leading order (NLO) corrections. In addition, the LO BFKL calculations have 
restricted phenomenological applications because, e.g., the running of the QCD coupling 
constant as is not included, and the kinematic range of validity of LO BFKL is not known. 

Recently the NLO corrections to the BFKL resummation of energy logarithms were 
calculated; see Refs. ||, |J] and references therein. The NLO corrections || |[ to the 
highest eigenvalue of the BFKL equation turn out to be negative and even larger than the 
LO contribution for as > 0.157. In such circumstances the phenomenological significance 
of the NLO BFKL calculations seems to be rather obscure. 

However, one should stress that the NLO calculations, as any finite-order perturba- 
tive results, contain both renormalization scheme and renormalization scale ambiguities. 
The NLO BFKL calculations [0, || were performed by employing the modified minimal 
subtraction scheme (MS) || to regulate the ultraviolet divergences with arbitrary scale 
setting. 

In this work we consider the NLO BFKL resummation of energy logarithms || |4[] 
in physical renormalization schemes in order to study the renormalization scheme de- 
pendence. To resolve the renormalization scale ambiguity we utilize Brodsky-Lepage- 
Mackenzie (BLM) optimal scale setting ||. We show that the reliability of QCD predic- 
tions for the intercept of the BFKL Pomeron at NLO when evaluated using BLM scale 
setting within non-Abelian physical schemes, such as the momentum space subtraction 
(MOM) scheme |7|, or the T-scheme based on T -> ggg decay, is significantly im- 
proved compared to the MS-scheme. This provides a basis for applications of NLO BFKL 
resummation to high-energy phenomenology. 

We begin with the representation of the MS-result of NLO BFKL ||, || in physical 
renormalization schemes. Although the MS-scheme is somewhat artificial and lacks a 
clear physical picture, it can serve as a convenient intermediate renormalization scheme. 
The eigenvalue of the NLO BFKL equation at transferred momentum squared t = in 
the MS-scheme || |J can be represented as the action of the NLO BFKL kernel (averaged 
over azimuthal angle) on the LO eigenfunctions [Q\j 1 Q\)~ l l 2+lu ||: 
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where 

X l{v) = 2^(1) - V(V2 + - ^(1/2 - w) 

is the function related with the LO eigenvalue, ip = T'/T denotes the Euler ^-function, 
the z/- variable is conformal weight parameter ||, Nq is the number of colors, and Qi^ are 
the virtualities of the reggeized gluons. 

The calculations of Refs. |3], |J allow us to decompose the NLO coefficient r 
(P) into /3-dependent and the conformal (/5-independent) parts: 



where 



and 



MS 



(2) 



' 4 



15 
2 XM ~ 3 



(3) 



con f I \ 
MS V I 



N, 



c 



2_ 4 



+ 



7T 



7r 2 sinh(7rz/) 
2^cosh 2 (7rz/) 

7T 3 



1 + 



iV>\ ll + 12z/ 2 



COsh(7T^) 



16(1 + z/ 2 ) 
-6C(3)+4^) 



with 



V?(z/) 



1 cos(z/ln(x)) 

Cl (1 + X)y/X 



~tt ~ Li 2 (x) 
o 



Li 2 (a;) 



eft 



ln(l - t) 



(4) 



(5) 



Here /3 = (ll/3)iV c - (2/3) JV> is the leading coefficient of the QCD /^-function, N F 
is the number of flavors, ((n) stands for the Riemann zeta- function, Li 2 (x) is the Euler 
dilogarithm (Spence-f unction). In Eq. (^) Np denotes flavor number of the Abelian 
part of the gg — ► qq process contribution. The Abelian part is not associated with the 



running of the coupling [ID] and is consistent with the correspondent QED result for the 
e + e~ cross section 11 . 

which is related to the running of the cou- 



The /^-dependent NLO coefficient r jfg{^) 



pling, receives contributions from the gluon reggeization diagrams, from the virtual part 
of the one-gluon emission, from the real two-gluon emission, and from the non-Abelian 



part [10j of the gg — > qq process. There is an omitted term in r^-(z/) proportional to 



MS K 



X^iy) which originates from the asymmetric treatment of Q% and Q2 and which can be 
removed by the redefinition of the LO eigenfunctions ||. 

The NLO BFKL Pomeron intercept then reads for Nq = 3: || 
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Physical renormalization schemes provide small and physically meaningful pertur- 
bative coefficients by incorporating large corrections into the definition of the coupling 
constant. One of the most popular physical schemes is MOM-scheme 0, |J, based on 
renormalization of the triple-gluon vertex at some symmetric off-shell momentum. How- 
ever, in the MOM-scheme the coupling constant is gauge- dependent already in the LO, 
and there are rather cumbersome technical difficulties. These difficulties can be avoided 
by performing calculations in the intermediate MS-scheme, and then by making the tran- 
sition to some physical scheme by a finite renormalization J/J. In order to eliminate the 
dependence on gauge choice and other theoretical conventions, one can consider renor- 
malization schemes based on physical processes ||, e.g., V-scheme based on heavy quark 
potential. Alternatively, one can introduce a physical scheme based on T — ► ggg decay 
using the NLO calculations of Ref. 
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A finite renormalization due to the change of scheme can be accomplished by a trans- 
formation of the QCD coupling 



a s -> a s 
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where T is some function of Nc, Np, and for the MOM-scheme, of a gauge parameter £. 
Then the NLO BFKL eigenvalue in the MOM-scheme can be represented as follows 
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For the transition from the MS-scheme to the MOM-scheme the corresponding T-function 
has the following form |7j]: 
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where I = -2 J* dx \n(x)/[x 2 - x + 1] ~ 2.3439. 
Analogously, one can obtain for the V-scheme 
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and by the use of the results of Ref. |C| for the T-scheme: 

6.47 Ar 2.77^ 
T T = ———Nc - -— A). 



(12) 
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Scheme 
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r(0) = r co "/(0) + r^(0) 


r(0) 
(N F = 4) 
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7.471 - 1.281/3b 


-12.64 - 0.1020iV F - 1.214/3q 


-22.76 





^ = 1 


8.247 - 1.281/3b 


-11.87 - 0.1020iV F - 1.214/3b 


-21.99 


M 


^ = 3 


8.790 - 1.281/3b 


-11.33 - 0.1020iV> - 1.214^ 


-21.44 
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2 - 0.4167/5o 


-18.12 - 0.1020iV> - 0.3497/3o 


-21.44 
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6.47 - 0.923/? 


-13.6 - 0.102iV F - 0.856^ 


-21.7 



Table 1: Scheme-transition function and the NLO BFKL coefficient in physical schemes. 

One can see from Table [1] that the problem of a large NLO BFKL coefficient remains. 
The large size of the perturbative corrections leads to significant renormalization scale 
ambiguity. 

The renormalization scale ambiguity problem can be resolved if one can optimize the 
choice of scales and renormalization schemes according to some sensible criteria. In the 
BLM optimal scale setting ||, the renormalization scales are chosen such that all vacuum 
polarization effects from the QCD /3-function are resummed into the running couplings. 
The coefficients of the perturbative series are thus identical to the perturbative coefficients 
of the corresponding conformally invariant theory with (3 = 0. The BLM approach has 
the important advantage of resumming the large and strongly divergent terms in the 
perturbative QCD series which grow as n![as/3o] n , i-e., the infrared renormalons associated 
with coupling constant renormalization. The renormalization scales in the BLM approach 
are physical in the sense that they reflect the mean virtuality of the gluon propagators 

BLM scale setting |J can be applied within any appropriate physical scheme. In 
the present case one can show that within the V-scheme (or the MS-scheme) the BLM 
procedure does not change significantly the value of the NLO coefficient r(u). This can 
be understood since the V-scheme, as well as MS-scheme, are adjusted primarily to the 
case when in the LO there are dominant QED (Abelian) type contributions, whereas in 
the BFKL case there are important LO gluon-gluon (non-Abelian) interactions. 

Therefore, from the point of view of BLM scale setting, one can separate QCD pro- 
cesses into two classes specifying whether gluons are involved to the LO or not. Thus 
one can expect that in the BFKL case it is appropriate to use a physical scheme which is 
adjusted for non-Abelian interactions in the LO. One can choose the MOM-scheme based 
on the symmetric triple-gluon vertex |7], [| or the T-scheme based on T — > ggg decay. The 
importance of taking into account this circumstance for vacuum polarization effects one 
can be seen from the "incorrect" sign of the /?o-term for rjjg in the unphysical MS-scheme 
(Eq- ©)■ 
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Scheme 


(JVf = 4) 


-l=u BLM (Q*,0) 


= 1 GeV 2 


Q 2 = 15 GeV 2 


Q 2 = 100 GeV 2 


M 

M 


e = o 


-13.05 


0.134 


0.155 


0.157 


e = i 


-12.28 


0.152 


0.167 


0.166 


^ = 3 


-11.74 


0.165 


0.175 


0.173 


T 


-14.01 


0.133 


0.146 


0.146 



Table 2: The NLO BFKL Pomeron intercept in the BLM scale setting within non-Abelian 
physical schemes. 



Adopting BLM scale setting, the NLO BFKL eigenvalue in the MOM-scheme is 
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The /^-dependent part of the tmom{v) defines the corresponding BLM optimal scale 
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Taking into account the fact that Xl{v) 
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Q 2 (4exp[2(l + 2J/3) - 5/3]) ~ Q 2 127. Note that 



Qblm 2 (v) contains a large factor, exp[-4T^ OM /^ ] = exp [2(1 + 27/3)] ~ 168, which 



reflects a large kinematic difference between MOM- and MS- schemes ||T3|, g|, even in an 
Abelian theory. 

Analogously one can implement the BLM scale setting in the T-scheme (Table 0). 

Figs. |l] and give the results for the eigenvalue of the NLO BFKL kernel. We have used 
the QCD parameter A = 0.1 GeV which corresponds to as = 47r/[/? ln(Q 2 /A 2 )] ~ 0.2 at 
Q 2 = 15 GeV 2 . Also, the generalization [|l~4| , |I~5| of the (3 -function in the running coupling 
and of flavor number for continuous treatment of quark thresholds has been used. 

One can see from Fig. [L], that the maximum which occurs at non-zero v is not as 
pronounced in the BLM approach compared to the MS-scheme, and thus it will not serve 
as a good saddle point at high energies. 

One of the striking features of this analysis is that the NLO value for the intercept 
of the BFKL Pomeron, improved by the BLM procedure, has a very weak dependence 
on the gluon virtuality Q 2 . This agrees with the conventional Regge-theory where one 
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Figure 1: //-dependence of the NLO BFKL eigenvalue at Q 2 = 15 GeV 2 : BLM (in MOM- 
scheme) - solid, MOM-scheme (Yennie gauge: £ = 3) - dashed, MS-scheme - dotted. LO 
BFKL (a s = 0.2) - dash-dotted. 



expects an universal intercept of the Pomeron without any Q 2 - dependence. The minor Q 2 - 
dependence obtained, on one side, provides near insensitivity of the results to the precise 
value of A, and, on the other side, leads to approximate scale and conformal invariance. 
Thus one may use conformal symmetry || [16| for the continuation of the present results 
to the case t ^ 0. 

Therefore, by the applying of the BLM scale setting within the non-Abelian physical 
schemes (MOM- and T- schemes) we do not face the serious problems [17, 18, |19[ which 
were present in the MS-scheme, e.g., oscillatory cross section disbehavior based on the 
saddle point approximation [ 17 1 , and the somewhat puzzling analytic structure [ 18 1 of the 
MS-scheme result [§, g. 

Now we will briefly consider NLO BFKL within other approaches to the optimization 
of perturbative theory, namely, fast apparent convergence (FAC) and the principle of 
minimal sensitivity (PMS) . 

By the use of the FAC |2(| one can obtain 



ufac(Q 2 , v) = N c Xl{v) 



as(Q 2 FAc(")) 



Ti 
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Figure 2: Q 2 -dependence of the BFKL Pomeron intercept in the NLO. The notation is as 
in Fig. 1. 



Q 



FAC\ 



Q 2 exp 



A 



-r(u) 



In the MS-scheme at v = 0, = 0.33 - 0.26 for Q 2 = 1 - 100 GeV 2 . However, 

the NLO coefficient r(u), and hence, the FAC effective scale, each have a singularity at 
vq ~ 0.6375 due to a zero of the xl (^-function. 

In the PMS approach |2l| the NLO BFKL eigenvalue reads as follows 
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where the PMS effective coupling olpms is a solution of the following transcendental 
equation 
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with C = /V( 4 A)) and A = 102 - 38iV>/3. At u = one obtains in the MS-scheme 
upms = 0.23 — 0.20 for Q 2 = 1 — 100 GeV 2 , but, by the same reason as in the FAC case, 
the PMS effective coupling has a singularity at v^. Thus, the application of the FAC 
and PMS scale setting approaches to the BFKL eigenvalue problem lead to difficulties 
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with the conformal weight dependence, an essential ingredient of BFKL calculations. The 
unphysical behavior of the FAC and PMS effective scales for jet production processes has 
been noted in Refs. 



22 



Before making conclusions a few remarks are in order. 

(i) Since the BFKL equation can be interpreted as the "quantization" of a renormal- 
ization group equation (TB| , it follows that the effective scale should depend on the BFKL 
eigenvalue u, associated with the Lorentz spin, rather than on v. Thus, strictly speaking, 
one can use the above effective scales as function of v only in "quasi-classical" approxi- 
mation at large-Q 2 . However, the present remarkable Q 2 -stabilty leads us to expect that 
the results obtained with LO eigenf unctions may not change considerably for t 7^ due 
to the approximate conformal invariance. This issue will be discussed in more detail in 
the extended version of this work [^3[ . 



(ii) There have been a number of recent papers which analyze the NLO BFKL pre- 
dictions in terms of rapidity correlations [24]], t-channel unitarity p5| , angle-ordering f2~6fl , 
double transverse momentum logarithms [27J and BLM scale setting for deep inelastic 
structure functions 



2N 



A discussion of these topics within our approach will be deferred 
toRef. 0. 

To conclude, we have shown that the NLO corrections to the BFKL equation for the 
QCD Pomeron become controllable and meaningful provided one uses physical renor- 
malization scales and schemes relevant to non-Abelian gauge theory. BLM optimal scale 
setting automatically sets the appropriate physical renormalization scale by absorbing the 
non-conformal /^-dependent coefficients. The strong renormalization scheme dependence 
of the NLO corrections to BFKL resummation then largely disappears. This is in contrast 
to the unstable NLO results obtained in the conventional MS-scheme with arbitrary choice 
of renormalization scale. A striking feature of the NLO BFKL Pomeron intercept in the 
BLM approach is its very weak Q 2 -dependence, which provides approximate conformal 
invariance. The new results presented here open new windows for applications of NLO 
BFKL resummation to high-energy phenomenology. 
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